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Abstract 

In the first part of the paper we present the dressing method which generates multi-soliton solutions 
to integrable systems of nonlinear partial differential equations. We compare the approach of Neuge- 
bauer with that of Zakharov, Shabat and Mikhailov. In both cases we discuss the group reductions 
and reductions defined by some multilinear constraints on matrices of the linear problem. The second 
part of the paper describes the soliton surfaces approach. The so called Sym-Tafel formula simplifies 
the explicit reconstruction of the surface from the knowledge of its fundamental forms, unifies various 
integrable nonlinearities and enables one to apply powerful methods of the theory of solitons to geomet- 
rical problems. The Darboux-Bianchi-Backlund transformation (i.e., the dressing method on the level 
of soliton surfaces) reconstructs explicitly many classical transformations of XIX century. We present 
examples of interesting classes of surfaces obtained from spectral problems. In particular, we consider 
spectral problems in Clifford algebras associated with orthogonal coordinates. Finally, compact formulas 
£Sj ■ for multi-soliton surfaces are discussed and applied for the Localized Induction Equation with axial flow. 

^t- . 

■ 1 Introduction 

rn 

An integrable system can be represented as intcgrability conditions for a linear problem (a system of linear 
partial differential equations containing the spectral parameter). In this paper we confine ourselves to linear 
problems of the form 

= U k * , (1) 

where x 1 ,x 2 are independent variables, := d&/dx k and U\, U 2 are complex nxn matrices which depend 
in the prescribed way on the dependent variables (shortly, "soliton fields" ) , also through their derivatives or 
integrals, and on the so called spectral parameter A. In practice Uk are meromorphic (usually rational) 
functions of A. 

Considering the system (1) as equations for an unknown function $ one can easily see that a non-trivial 
solution exists if and only if the following intcgrability condition holds: 

U 1 , 2 -U 2 , 1 +[U 1 ,U 2 }=0 . (2) 

The "zero curvature condition" (2), considered as an identity with repect to A, is equivalent to a system 
of nonlinear partial differential equations for soliton fields. The system (2) has many interesting properties 
which justify to call it integrable [2, 48, 59, 85]. One of the properites, the existence of the Darboux- 
Backhmd transformation, is discussed thoroughly in this paper. 

The system (1) has n independent vector solutions. It is convenient to consider them as columns of a 
non-degenerate nxn matrix called the fundamental solution of (1). In the sequel we always assume that <f>, 
sometimes called the wave function, is the matrix. 

Remark 1 A linear combination of vector solutions is a solution as well. In other words, if 4> and $' are 
fundamental solutions of (1) then there exists a constant non- degenerate matrix C such that $' = $C 
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As an example consider the nonlinear Schrodinger equation: 



iq,2+fq, 1 l+2q\q\ 2 =0 , (3) 
where q — q(x 1 ,x 2 ) € C is the soliton field. The associated linear problem reads [2, 85], 

$ 7l = (i\a 3 + Q)$ , $, 2 = (-2U 2 (Tg -2\Q + R)<P , (4) 



where 03 := 




The Darboux-Bianchi-Backlund (DBB) transformation is a gauge transformation which preserves 
the form of the linear problem. On the level of soliton fields this transformation simply adds a soliton 
solution. Usually the name of Bianchi is not mentioned in this context. However, because of his enormous 
contribution to this field it seems justified to name after him at least the corresponding transformation on 
the level of soliton surfaces (see Section 4). 

There are several approaches to the construction of multi-soliton solutions by gauge transformations. Let 
us mention the Zakharov-Shabat dressing method [51, 85, 89], the Neugebauer approach [61, 65], Darboux 
transformation for linear differential operators [58, 59], the method of Its based on axiomatics of the wave 
functions [43] and the approach of Gu [40] . The methods have been discovered to some extent independently. 
The fundamental concept of the "dressing" of linear operators proposed by Zakharov and Shabat [88, 89] was 
probably the starting point of all these researches. The details however are different and very often authors 
do not seem to be aware of the results of the other groups. In this paper we focus on the Neugebauer and 
Zakharov-Shabat methods. 

The dressing method of Zakharov and Shabat [85, 86, 88, 89] is a general scheme to solve integrable 
nonlinear equations. Consider a gauge transformation 

$ = (5) 

defined by a matrix D = D(x 1 ,x 2 ; A). Then, obviously, $ satisfies a linear system (1) with matrices Uk 

Uk^D^D^+DUkD- 1 , fc = l,2. (6) 

If we are able to construct a matrix D such that Uk given by (6) are of exactly the same form as Uk then D 
is called the Darboux matrix [18, 52, 59]. In other words, the dependence of Uk on A and on soliton fields 
should be exactly the same as that of Uk (of course, the soliton fields entering Uk are in general different 
than those entering Uk)- The existence of the Darboux matrix is one of the criterions of intcgrability. 

It turns out that the structure of matrices Uk can be characterized completely in terms of some simple 
properties. The most important information is contained in the singularities of Uk- The matrices Uk are 
assumed to be meromorphic in A and their poles are given. 

In this paper we show that the construction of the Darboux matrix is practically algorithmic (see also 
[18]). The crucial point is to notice all relevant algebraic and group properties of the associated linear problem 
(especially to identify the reduction group) and then to apply appropriate general theorems. In particular, 
we present multilinear constraints which are usually invariant with respect to DBB transformation. Taking 
them into account we can avoid some cumbersome calculations and our construction assumes a more elegant 
form. 

It is reasonable to consider as equivalent the Darboux matrices which generate the same DBB trans- 
formation on the level of soliton fields. In other words, D\ and D2 are equivalent iff the corcsponding 
transformations (6) are identical. 

Remark 2 The Darboux matrix D is equivalent to fD where f is any scalar function independent on x , x , 
i-e.f = f(X). 

The methods of constructing the Darboux matrix are based on various approaches to matrices with 
rational coefficients. The following example shows that problems of that kind can be interesting even in the 
well known scalar case. Namely, we proceed to the decomposition into partial fractions. The procedure (used 
in the integration of rational functions) is standard and algorithmic but usually is considered as cumbersome 
and boring. However, using local parameters in the neighbourhood of singularities one obtains a very effective 
method. 
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Let us show the main idea of the approach on the example of the function F(x) = (2x + l) w x 3 . Taking 



The crucial point is that h(x) = 0. Indeed, by construction h is bounded and regular both in x — and 
x = —1/2. Hence, by the Liouville theorem we have h(x) = const. The constant is zero by the obvious 
boundary condition: h — > for \x\ — > oo. The presented approach simplifies computations in a striking way. 
We can integrate very quickly any rational function with known poles. 

The Liouville theorem is the corner stone of the dressing method (as pointed out by Its [43]): local 
behaviour of analytic functions (e.g., around singularities) can define them globally. However, the matrix 
case is much more complicated (and more interesting) than the scalar one (compare [37, 68]). 

In this paper we apply methods of the theory of integrable systems to the geometry of surfaces immersed 
in Euclidean spaces. The main idea of the soliton surfaces approach consists in associating with a given 
integrable system of nonlinear partial differential equations a class of surfaces (or manifolds) immersed in the 
Lie algebra of the corresponding linear problem using the so called Sym-Tafel formula (20) (see [79], compare 
also [21, 34, 42]). In this way we can unify in a natural way various nonlinear models (spins, vortices, chiral 
fields, etc.). 

We present several examples, mostly immersions in R 3 : surfaces swept out by various vortex motions, 
surfaces of constant Gaussian curvature, constant mean curvature surfaces, Bianchi surfaces, isothermic 
surfaces. The Sym-Tafel formula can be also applied to manifolds immersed in spaces of higher dimensions 
like orthogonal coordinates in R™ [23] and n-dim. space forms immersed in R 2 ™ -1 [22]. 

The Darboux-Bianchi transformations (DBB transformations on the level of surfaces) coincides, as 
a rule, with classical transformations studied by Luigi Bianchi and other great geometers of XIX century 



There are several approaches to the construction of the Darboux matrix. All these methods are closely 
related but usually not much attention is given to this fact. The discussion of the relationship between the 
Matveev method and the Zakharov-Shabat approach is a positive exception ([6], see also [59]). In this paper 
we compare the Zakharov-Shabat method with the approach of Neugebauer. 

The Neugebauer-Meinel approach 

An especially effective method to construct iV-fold Darboux-Backlund transformation has been proposed by 
Neugebauer ([61, 65, 66], see also [40]). The Darboux matrix is assumed to be a polynom in A, 



(polynomial Darboux matrices will be denoted by T). Any matrix D rational in A can be replaced by an 
equivalent matrix T of the form (7) (compare Remark 2 where for / we should take the least common 
denominator). One can consider other equivalent matrices, e.g., polynoms in 1/A. 

From the elementary linear algebra we know that T _1 = (det T)~ 1 T', where T" is the matrix of cofactors 
of T. Obviously T" is a polynom in A. Therefore, if Uk are rational functions of A, then Uk given by (6) are 
rational as well. What is more, the only candidates for poles of Uk are poles of Uk and zeros of detT. 

Lemma 3 // TiUk = then the determinant of the Darboux matrix does not depend on x 1 , x 2 . In particular, 
the zeros of detT are constant. 

The proof is based on a very useful identity, (log det A), x — Tr(A, x A -1 ), which holds for any nondegener- 
ate matrix function A = A(x). The conditions TrJ/fc = 0, assumed throughout the paper, are not particularly 
restrictive: practically all linear problems associated with integrable systems are tracelcss. 

The necessary condition for T to be a Darboux matrix is the requirement that Uk have no poles in 
zeros of detT. We assume that det Co 7^ and all zeros of detT (denoted by Ai) are simple and pairwise 




[8, 49, 79]. 



2 The Darboux-Backlund transformation 




(7) 



3=0 
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different. The total number of zeros is Nn. We assume also that Ui, U2 and $ are regular in A*. Because 
Uk = $,fe = $,fe $'(det $ dct T) -1 , then the condition for U k to have no pole in A = \ is given by 

$,k(X i )$'(X i )=0 . (8) 

We proceed to construct $ in order to satisfy the equations (8). The matrix $(Aj) is degenerated. Hence 
there exists a vector pi such that $(Aj)p, = 0. Then, by (1) it follows that (Xi)pi = 0. Differentiating 
the equation $(Aj)pj =0we obtain $(Aj)pi,fc = 0. Therefore, taking into account that zeros Aj are simple, 
we obtain that pi,k is proportional to pi. In fact, using an additional freedom (pi are defined up to a scalar 
factor) we can choose pi to be constant vectors. 

Corollary 4 If Aj is a simple zero ofdetT, then there exists a constant vector pi 7^ such that $(Aj)p, = 
$,fc (Ai)pi = 0. 

The corollary implies the equation (8). Indeed, one has only take into account the following fact of the 
elementary linear algebra (we suggest to prove it as an exercise) [61]. 

Lemma 5 If there exists a vector p 7^ such that Xp = and Yp = (X,Y are matrices) then YX' = 0. 

Finally, we can treat constants X t e C and Pi <G C" (i = 1, . . . , nN) as prescribed parameters. Then the 
equations 4>(Aj)pj = will determine the coefficients Cj of the Darboux matrix (7). Namely, the system 

T(A fe )$(A fe )p fe = , {k = l,...,Nn) , (9) 

is equivalent to Nn 2 scalar equations for (N + l)n 2 scalar coefficients. The remaining freedom corresponds 
to the gauge freedom. For instance we can treat C as still undetermined. Then matrices C\,. . . ,Cjv are 
uniquely determined provided that the complex parameters Aj. and constant complex vectors pk are given. 

The degenerate case (multiple zeros of detT) is more complicated. If there exist several vectors 
satisfying the equation $(\j)pj = 0, then they also can be chosen to be constant. 

To obtain explicit solutions of a given nonlinear system by the presented method we have to know 
explicitly at least one solution of the linear problem (1) . Then we can obtain in an algebraic way a sequence of 
other solutions which can be interpreted as a "nonlinear superposition" of the given ( "background" ) solution 
and some number of solitons. In particular, applying the transformation (7) to the trivial background we 
obtain "pure" iV-soliton solutions. 



The Zakharov-Shabat approach 

In the papers of Zakharov, Shabat and Mikhailov another representation of the Darboux matrix was used 
[63, 85, 86, 87], 

where I is the unit matrix, N,Ak,Bk depend on x 1 ,x 2 and \k,Pk are constants (assumed to be pairwise 
different). 

The similar form of D and D^ 1 is a restriction on Ak and Bk- The condition DD^ 1 = I is equivalent to 

= 0, (j + V^V ) B k =0, (11) 
(k = 1, . . . , N). In the case N = 1 ("1-soliton case") this system can be easily solved to give 

D = "(' + feH- °-'-(, + !£±ry>, (12) 

where P is a projector (i.e. P 2 = P). 

Theorem 6 (Zakharov, Shabat [89]) If P is given by the formulas: 

kerP = ^(X^Vker , imP - ${m)V im , (13) 

where Vk er and Vi m are constant vector spaces such that Vk er © Vim — C", then the transformation (6) with 
D given by (12) preserves the divisors of poles of matrices Uk- 
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Every projector is completely characterized by its image and kernel. The matrix of the projector P is 
explicitly given by 

P = {0 1 ^i)V m }{<i>(X 1 )V ker ,Mp 1 )V m }- 1 , 

where in the brackets we have nxn matrices. We use the same notation to designate a subspace of R" and 
a matrix representing it (the columns of the matrix span the considered vector space) . 

Multiplying D given by (10) by the common denominator we obtain an equivalent polynomial matrix 
like (7). However (10) and (7) are not strictly equivalent. One should remember about restrictions (11). An 
equivalent matrix D of the form (10) exists only for some particular polynomial matrices T. 

3 Algebraic representation of the linear problem 

The techniques presented above guarantee that a linear problem subject to the Darboux-Backlund trans- 
formation is transformed into a linear problem of the same analytical dependence on A. However, a given 
integrable system usually needs much more restrictions imposed on the associated linear problem. The most 
important reductions consist in confining Uk to some Lie algebra. Then $ is confined to the corresponding 
group and, obviously, D is confined to the same group. An example is in order. Let l^(A) = -U k {\) 
for k = 1,2 (this property holds in the case of the linear problem (4) where the coefficients by powers 
of A arc su(2)-valucd). Then one can check that ((&(\))- 1 ), k = U k (X)(& (X))- 1 and, by virtue of Re- 
mark 1, ($t(A))-! = $(A)C(A) where C(A) does not depend on x 1 ,x 2 . Elementary calculations show that 
C^(A) = C(A). In fact C depends only on the initial condition $(xJ,Xg;A) and by an appropriate choice 
of the condition we can put C = I. Thus <I> satisfies $ 1 '(A)(f>(A) = /. In general, if U k are confined to 
some loop algebra then, by an appropriate choice of initial conditions, $, $ and D can be restricted to the 
corresponding group. 

The reduction groups had already been introduced in the pioneering paper of Zakharov and Shabat 
([89], see also [85]). Then the subject was developed by Mikhailov [62, 63, 86, 87]. He classified a large class 
of reductions considering groups of 4 types [62]. 

The restrictions on the matrix D given by (12) has been discussed in detail in [18]. Actually, because of 
geometrical applications, it is convenient to consider the unimodular Darboux matrix V [18, 79] (by Lemma 3 
it is clear that in the case of traceless linear problems one can always impose the condition det£> = 1). To 
confine the Darboux matrix (12) to the group SL(n, C) it is sufficient to multiply it by the scalar coefficient 
/ = (dctD) 1 /"; 

where detA/ r =l, d=dim(imP) and Ai ^ fi\. The matrices (12) and (14) are equivalent (compare Remark 2). 
The paper [18] contains the full list of restrictions on on Ai, /xi, /, P,Af implied by the Mikhailov reductions 
provided that det M ^ 0. 

Here we are going to present briefly some results concerning Darboux matrices in the polynomial form 

(7). 

Let Uk(— A) = JU(X)J 1 where J can depend on A. One can easily see that J has to satisfy J(— A) J(A) = 
0(X)I where 9 is a scalar function. The wave function and the Darboux matrix can be confined to the group 
$(-A) = J$(A)J-\ T(-A) = JT(A)J" 1 . Therefore dctT(-A) = detT(A) which means that zeros of 
detT(A) appear in pairs Xj = — X k . The corresponding constant eigenvectors pj and p k satisfy (fr(Xk)pk = 
and $(Aj)pj = which implies $(Afe) J(— X k )pj = 0. If the zeros Aj are simple then pj = J(X k )pk- 

The unitary reductions Ul(X) = —HUk(X)H~ 1 are more complicated. If we admit a (rational) dependence 
of H on A then H^(X) = rj(X)H(X) where r\ is a scalar function. One can easily check that iJ($'l'(A)) _1 satisfies 
the system (1) provided that $ is a solution. <E> is subject to an analogical constraint. Thus we can derive 
the condition 

T(X) = k(X)H (THX))- 1 H- 1 , (15) 

(k is a scalar function of A) which is sufficient for T to be the Darboux matrix. We point out that now the 
simplest choice k(X) = 1 is not possible because it is incompatible with the form (7). If H is rational then k 
has to be rational as well. From (15) we obtain 

(fc(A))™ = det T(X) det T(A) . (16) 
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Therefore k(X) is a polynom with real coefficients. We confine ourselves to the case of k(X) without real 
roots, i.e. 

k(X) = (A - Ai)(A - AO ... (A - Aat)(A - X N ) 
where N is the degree of the polynom T(A). Then 

detT= (A- A 1 )"- dl (A-A 1 ) dl ...(A- A w )™- d "(A-A Ar ) djv . 

where dk are some integers. Note that for n > 2 the zeros of det T are, as a rule, degenerated. If Xk are 
pairwise different, then dividing T by (A— Ai) . . . (A— Ajy) we obtain the matrix D (equivalent to T) bounded 
for A — > oo: 



D _ T(X) _ ? ■ y S * 



(A-A 1 )...(A-A A r) u ^ X-X k 

where Sk are some matrices dependent on x 1 ,^ 2 . From (16) we can deduce that the matrix So is not 
degenerated. Therefore D is exactly of the form (10). Using (15) we compute D~ x and, as a result we have 
lik = Afe. The case of k(X) with real roots is left as an exercise. We just mention that this case is possible 
only for even n. 

General conclusion is very interesting: in the important case of unitary reductions the Darboux matrix 
can always be represented in the form (10). 



Multilinear invariants of the Darboux matrix 

Let Ao be a pole of matrices U\ and U2 of the order J and K respectively. The matrices U := (A — Xq) j Ui 
and W := (A — Xq) k U2 are holomorphic in the neighbourhood of A : 

00 00 

= (A -A )\ W = Y,w i {X-Xo) i . 

1 = 4=0 

We assume that Uj and Wj are sl(n, C)-valued functions of x , x 2 . In many interesting cases including the 
equation (3) the only singularity is Ao = 00, i.e. U\ and U2 are polynomials in A while U, W are polynomials 
in 1/A. 

The transformation (6) for U, W reads: 

U = (A - A ) J L>,i D- 1 + DUD- 1 , W = (A - X Q ) K D, 2 D' 1 + DWD' 1 . (17) 

Considering the leading terms in equations (17) we can obtain a number of invariants of DBB transformation. 
In the sequel we assume J < K and D is regular at A = Ao. 

Linear invariants. Consider aU + j3W where a — a(x 1 ,x 2 ) and (3 — ^(x 1 ,^ 2 ) are given functions. We 
have 

aU + pW = a(X - A ) J L>,i D^ 1 + /3(A - X Q ) K D, 2 D' 1 + D(aU + (3W)D- 1 . 

Let L < J < K. If aU + (3W has zero of L-th order in A = Ao then aU + f3W has zero of at least the same 
order. Thus we have shown the following result [18]. 

Proposition 7 If X = const then any system of L + 1 linear constraints (L < J) of the form 

au m + pw m = , (m = 0, 1, . . . , L) , 
where a — a(x 1 ,x 2 ), (3 — ^(x 1 ,^ 2 ) are given functions, is preserved by DBB transformation. 

The proposition can be extended for L = J provided that the Darboux matrix satisfies the condition 
D,i (Ao) = (in the case Ao = 00 it means that the normalization matrix does not depend on x 1 ). 

In the case of the nonlinear Schrodinger equation we have (by inspection): u + 2w — and u\ + 2w\ = 
(see (4)). These constraints are invariants of DBB transformation provided that N = N{x 2 ). 

Bilinear invariants. Let the center dot denotes an invariant scalar product in sl(n, C), namely a ■ b := 
Tr(a&). Consider the transformation of W ■ W: 

W ■ W = W ■ W + 2(A - Ao) x (D _1 A2 ) • W + (A - X^iD^D^ ) ■ (D~ 1 D, 2 ) . 
On the right hand side the leading terms (of order less than K) are contained in W ■ W. 
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Proposition 8 If X — const then the constraint 

m 

^2 W k ■ Wm-k = 7m (0<m<K-l) 
k=0 

(where j m — j m (x 1 7 x 2 ) is a given function) is preserved by DBB transformation. 

The proposition holds for m = K if and only if (D^ 1 (X )D,2 (Ao)) • w — (in the polynomial case 
(Ao = oo) this condition reduces to N ,1 = 0). 

The multilinear invariants can be used to obtain reductions of the linear problem (1). The reductions 
are fixed by the choice of a, (3 and 7 TO . 

Similar results can be obtained considering the expansions of U ■ U and U ■ W around A = Ao. One can 
expect to get analogical results for Tr(UU ...U), Tr(WW ...W), Tr(UW . . . W) etc. 

The linear problem as a system of algebraic constraints on two matrices 

Let us summarize the results of the preceding sections. The knowledge of the linear problem like (1) is 
crucial to solve a given integrable system. One can describe the linear problem parameterizing it explicitly 
by soliton fields, their derivatives and integrals. However, if the the construction of the Darboux matrix is 
concerned, much more convenient description is by a system of relatively simple (non-differential) constraints 
(algebraic representation of the linear problem [18, 19]): matrices XJ\ and Ui are rational functions of 
A (with prescribed poles) restricted to some loop algebra (in other words, their coefficients by powers of A arc 
restricted to a Lie algebra) and some multilinear constraints are imposed on the coefficients by powers of A. 
Considering nonisospectral linear problems [13] we should also prescribe an explicit dependence of variable 
spectral parameter on a constant parameter [18]. 

The algebraic representation of the linear problem allows us to introduce a very useful definition of 
the Darboux matrix. The matrix D is called the Darboux matrix for the linear problem (1) if the 
transformation (6) preserves the system of algebraic constraints equivalent to this linear problem [18]. 

It is very convenient to divide the construction of the Darboux matrix into two separate steps ([18, 19], 
similar idea can be found also in [43]). First, we represent the linear problem under consideration as a system 
of algebraic constraints on two matrices. This step may turn out to be rather non-standard (i.e. each case 
has to be treated individually) but involves no cumbersome calculations. Usually the algebraic constraints 
can be found by straightforward inspection of the given linear problem. The crucial point is to find the 
reduction group. Second, we derive the constraints on the Darboux matrix imposed by the requirement that 
the matrix has to preserve the algebraic constraints. The second step is technically much more difficult but 
there exist many general propositions which make the construction practically algorithmic (the 1-soliton case 
is discussed in detail in [18]). 

4 Soliton surfaces approach 

We proceed to apply the soliton theory to the geometry of surfaces immersed in Euclidean spaces. Let us 
consider a surface £ C R 3 equipped (locally) with coordinates x 1 , x 2 , and defined by the position vector 
r = r(x 1 ,x 2 ) e R 3 . The first fundamental form ("metric tensor") is given by I := dr ■ dr (the center 
dot means a scalar product in the ambient space R 3 ). Therefore, in local coordinates, / := gijdx l dx^ 
(we assume the Einstein convention, i.e., summation over repeating indices), where gij := r,i-r,j. The 
matrix inverse to (g^) will be denoted by (g^) and g := det(gij). Let n be the unit normal vector, 
n := g _1 / 2 r,ix r, 2 , where the cross denotes skew (or vector) product in R 3 . The second fundamental form 
is defined as II := —dn-dr, or equivalently, II — b i jdx l dx : > where bij := r,^ n. The matrix with coefficients 

'■= g lk bkj represents the "shape operator" [11]. The principal curvatures, k\ and k 2 , are defined as 
eigenvalues of {¥■). The corresponding eigenvectors ("principal directions") can also be characterized by the 
requirement to diagonalize both fundamental forms. The Gaussian curvature is defined as K := k\ki and the 
mean curvature as H := ^{ki + k 2 ). They can be expressed by the fundamental forms: K = det(bij) /det(gij), 
H = g^bij. In fact, K can be expressed solely in terms of gtj (the Gauss "Theorema Egregium"). 

Differentiating tangent vectors we obtain Gauss- Weingarten (GW) equations r,y = F*jr,fc +bijia , 
where are some x 1 , x 2 -dependent coefficients known as "Christoffel symbols". It is well known that they 
can be expressed entirely by g^, namely T^- = \g ks {gi S ,j +gsj,i —gij,s )• The GW equations can be rewritten 
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as kinematic equations of the orthogonal frame ei, e 2 , n, where e k form an orthogonal basis in the tangent 
space: 





fifc e 2 , (18) 



where Q\ and S^2 are so(3) matrices depending on gij and bij. It is convenient to use the well known 
isomorphism so(3) ~ su(2) (see, for instance, [28]) and to rewrite the equations (18) in terms of 2 x 2 
matrices. The explicit form of the matrices reads, for example, as follows (compare [21, 77]): 

( gnfc2fc-gi2&ifc l , I ,-p2 / a \ 
* r~ lk+ f lk \ 911 I* (* = i,2)- (19) 

Integrability conditions for GW equations can be derived from r,ij k = Y,n~j. They are expressed explicitly in 
terms of matrices £l k : 

and are equivalent to a system of 3 nonlinear equations for the coefficients g^, b^: the Gauss-Mainardi- 
Codazzi (GMC) equations. 



The Sym-Tafel formula 

Let matrices Ui, Ui for A e R assume values in a Lie algebra g (in other words, matrices U k are meromorphic 
functions of A with g- valued coefficients). Then one can confine $ to the corresponding Lie group which 
implies that the matrix r = r(x 1 ,x 2 ; A) given by 

r = (20) 

is g-valued. The formula (20) is known as the Sym formula or the Sym-Tafel formula (the final step in 
derivation of this formula was done by Tafel, see [72]). We assume that the Lie algebra g is equipped with 
a scalar product invariant with respect to adjoint transformations, (Ad<t,a\ Adq,b) — (a\ b) (in a matrix 
representation Ad$,a := <f> _1 a<f>). A well known example is the Killing-Cartan form, (a\ b) := Ty(ad a ° adb), 
which is non-degenerate for semi-simple algebras. In general, Tr(/(a)/(6)) defines an invariant bilinear form 
for any matrix representation / of the Lie algebra g, and can be used as a scalar product provided that 
the bilinear form is non-degenerate [39]. Then g can be identified with a pseudo-Euclidean space R m and 
the function r represents a A- family of parametric surfaces [83] ( "soliton surfaces" ) immersed in g = R m . 
Kinematics of r defines an integrable nonlinear model which can be interesting in itself [73, 79]. 

Starting from the formula (20) we compute tangent vectors r,k= &~ 1 Uk,\ $ and the metric tensor 
gij = (r,j| r,j ) = (Ui,\\ Uj,\ ) (the invariance of the scalar product was used). We obtained gij without solv- 
ing differential equations (1). It was sufficient to know only the matrices U\, U 2 - GW equations read r,/y = 
$ -1 (t/fc,Aj +[U k ,\ , Uj])& and GMC equations (r, ijk = r, ikj ) are given by $ _1 [J7 i)A , Uj, k -U k ,j +[Uj, U k ]}$ = 
0. 

The Sym-Tafel formula gives an interesting connection between the classical geometry of manifolds (with 
possible singularities) immersed in R m and the theory of solitons [79] . The soliton surfaces approach is very 
useful in construction of the so called "integrable geometries" [8, 15, 25, 79]. Indeed, any class of soliton 
surfaces given by (20) is integrable. Geometrical objects associated with soliton surfaces (tangent vectors, 
normal vectors, foliations by curves etc.) usually can be identified with solutions to some nonlinear models 
(spins, chiral models, strings, vortices etc.) [73, 79]. The transformation (5) can be immediately extended 
on the models associated with soliton surfaces. For instance, the transformation (5) applied to the position 
vector (20) assumes the form 

f = r + <p- 1 V- 1 V, x ^ (21) 
and we propose to call it the Darboux-Bianchi transformation. Usually r€jC sl(n, C) and we may 
consider the unimodular Darboux matrix (14). In the SU(n) case we have the additional restriction yUi = Ai. 
The transformation (21) for soliton surfaces immersed in su(n) is given by [73]: 

2ImAi ^ , / d \ ^ 

r = r+ [--P (22) 

I A- Ai| 2 \n J ' 

where Ai is a complex parameter and P is the hermitean projector (i.e., P 2 = P and P^ = P) onto the 
space &{\i)V im , where V im is a constant vector space of a given dimension d. The length of the segment 
f — r does not depend on x 1 , x 2 . 
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Geometry from spectral problems 

The Euclidean 3-dimensional space E 3 can be identified with the Lie algebra su(2) endowed with the scalar 
product given by the Killing-Cartan form 

a-b:=(a\b) = -2Tr(a6) . (23) 

Note that Tr(a 2 ) < for a e su(2). What is more, the skew product in R 3 can be identified with the 
commutator of su(2) matrices. The standard basis, orthonormal with respect to (23), is given by 

e i = -^(? n). e 2 = -iz(° "*], e 3 = -li(l ° ). (24) 



2 V 1 J ' z 2 \ i J ' ° 2 V -1 
In other words, wc identify 



1 



su(2)3-^l^ -y-«° ) <- (,, y ,,),E 3 . 

Below we present several examples of integrable geometries in E 3 . We start from given linear problems 
and use the formula (20). The normal vector to r is given by n = (— 2Tr([i7 1 , A , U 2 ,\ ] 2 )) _1 / 2 $ _1 [Ui,\ U 2 ,\}$. 
Then we compute the fundamental forms and the mean or Gaussian curvature of the resulting surface. 

In order to put the result into a more elegant form we have chosen a special dependence between A and 
the spectral parameter ( — ((X). In two cases (constant mean curvature surfaces and spherical surfaces) 
the linear problem depends on the spectral parameter through trigonometric functions sin/« and cosk, where 
k is proportional to A. It is equivalent to using the spectral parameter £ confined to the unit circle, i.e., 
sinK = Yi (C — an d cosk = \{C + ^)- Wc use notation x 1 = u, x 2 = v. 

• Pseudospherical surfaces 

$,i= (-Ce 3 - ip,ie 2 )$ , 
<P, 2 = (, 1 (costp e3 — simp ei)<P , 

where £ = exp(— RX). The compatibility conditions (2) are equivalent to the sine-Gordon equation 
¥?;12 = simp and the fundamental forms are given by 



(26) 



I = R 2 (( 2 du 2 + 2cosipdudv + (~ 2 dv 2 ) , 
II = 2R simp dudv . 

The Gaussian curvature is given by K = —R 2 (compare [73, 79]). 
Constant mean curvature surfaces 

$,i = (-(ef/ 2 + e -t/ 2 cos2/«)ei - \tp, 2 e 2 - er^l 2 sin2«e 3 ) $ , 
$, 2 = (-e-^/ 2 sin2« ei + ±<p tl e 2 - (e*>/ 2 - er^l 2 cos2/t)e 3 ) $ , 

where k := ^H~ 1 X. GMC eqs. are reduced to the single equation 

<p,n +<P,22 +4 sinhf^ = , (28) 



(27) 



known as the elliptic sinh-Gordon equation. The fundamental forms of (20) read 

/ = H- 2 e^(du 2 + dv 2 ) , 
II = H' 1 ({e* + cos2n)du 2 + 2 sm2n dudv + (e^ - cos2 K )dv 2 ) , 

which means that the surface has the constant mean curvature H [31]. 



(29) 
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• Spherical surfaces 

$,i = (—2 sinh^ cosk ei — \<p,i e 2 + 2 cosh^ sum 63) $ , 
$,2 = (—2 sinh^j sin/t ei + \ip,\ e 2 — 2 cosh^ cosk e 3 ) $ , 

where k := RX, and GMC eqs. are equivalent to (28). 

I = 2R 2 ((coshys + cos2k) du 2 + 2sin2ft dudv + (coshes — cos2k) dv 2 ) 
II = 2Rsmhip (du 2 + dv 2 ) . 

Therefore, r defines a spherical surface (K = R~ 2 > 0) [31]. 

• Bianchi surfaces 

*,i = (-C«e 3 - (p,i +^ sin^) e 2 ) $ , 
*,2= (j(bcosif e 3 -6sin</? ei) + £j±£sin<pe 2 ) $ , 
where p — f(u) + g(v) (/, g are given functions) and 

A fl-2Xg(v)\ 1/2 
C = ; , ow \ > ( A = const ) . 



(30) 



(31) 



(32) 



I + 2A/(u), 
is "variable spectral parameter" . We obtain 

I = (C,a ) 2 (a 2 dw 2 + 2a6C" 2 cosp dudv + b 2 C 4 dv 2 ) 

II = — 2(,\(~ 1 ab simp dudv . 



(33) 



The Gaussian curvature is K = — p 2 where p := 1+2 ^ M ) + i-2\g{v) • ^ or ^ = we have p = f + g 
and the formulas (33) assume the standard form: 

I = p 2 (a 2 du 2 + 2ab coslo dudv + b 2 dv 2 ) , 

V ; (34) 

II = 2pab sirup dudv . 

We recognize the fundamental forms of Bianchi surfaces in asymptotic coordinates (for more details 
see [18, 46, 53, 80]). 

Localized Induction Equations and multi-soliton surfaces in R 3 

Soliton surfaces in R 3 and several nonlinear models of physical importance [26, 55, 76, 78, 79, 27] are 
associated with su(2) algebra. The corresponding linear problems (known as su(2)-AKNS linear problems 
[2], see also [79]) are parameterized by analytic functions u> = ui(X). A typical example is given by the linear 
problem (4) of the nonlinear Schrodinger equation (oj(X) = — 2X 2 ). The function oj uniquely characterizes 
the asymptotic behaviour of multi-soliton solutions. The linear problem corresponding to the trivial solution 
(g = 0) has the form 

$0, x = iX<7 3 $ , <j> , t = iw(A)(T3$o , 

and can be solved easily 

$0 = cxp(iXx<7s + iu)(X)tai) . (35) 
The case w = — 2aX 2 — 4/3A 3 is associated with the Hirota equation 

iq, 2 +a(q,u +2| q \ 2 q) - i0(q,m +6| q | 2 g,i ) = , (36) 

where a, (3 are real constants and q = q(x 1 ,x 2 ) € C. The kinematics of the position vector to soliton surfaces 
(evaluated at A — 0) of the Hirota equation describes the motion of the single thin vortex filament in the so 
called localized induction approximation with axial flow [36, 44]: 

r, 2 = a(r, 1 xr, 11 )+/3(r, 111 +^(r, 11 ) 2 r, 1 ) , (r ;1 ) 2 = 1 , (37) 
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where r = r(x 1 ,x 2 ) G R 3 . The tangent vector S :— r,i solves the following spin model 



S, 2 = a(SxS, n ) + (3 (s, n 

The special case /3 = corresponds to Localized Induction Equation, nonlinear Schrodinger equation and 
continuum Heisenberg ferromagnt model respectively (see, for instance, [35, 50, 67, 73]). In the case a = 0, 
q G R the equation (36) is known as modified Korteweg-de Vries equation. Soliton surfaces are degenerated: 
the position vector r describes the evolution of a plane curve which has an interesting elastomechanical 
interpretation and admits so called "loop solitons" as special solutions [76] . 

In the su(2) case usually the formula (20) is used with with the factor i. Therefore, in this section 
we assume r = i<f> _1 $,A. Iterating the Darboux-Backlund transformation N times we obtain the following 
expression for TV-soliton surfaces (more precisely, 7V-soliton addition to any background surface) associated 
with su(2)-linear problems [14, 79]: 

A, ( 2ReS fc 2ImS fc j S fc j 2 - 1 \ 

r B +N = r B + 2^ d k l - T — ei - 2 e 2 + 2 e 3 , (38) 

fe=1 V I Sfe I + 1 I I + 1 l =k l + 1 / 

where the subscript B means "background" , are defined by (24) , dk are given by 

7 ImAfe 

d fc := 77 Tl2 - ( 39 ) 

| A - Afe | 

Afe are constant complex parameters, and Ek parameterize P k , orthogonal projectors onto 1-dim. subspaces 
of C 2 , defined by imPfe + i = $^^_ fe (A)$ B+k(^k)Pk where pk G C 2 are constant vectors. Namely, 

Pfe "i + |s fc K 2 i-|s| 2 

Trying to compute tb+tv one meets serious technical problems even for N — 2 (compare [55]). However, it 
is possible to simplify the problem and to express the formula (38) in terms of functions (k = 1, 2, . . . , N) 
defined by 

6, := ^ , where f ) := ^(^(A*) f f ) , (40) 

U k 2 \ u k2 ) \ Pk2 ) 

where Pki,Pk2 are components of pk (see [14]). Let us represent the complex function by 

£ k =■ cxp(Q fe - iotk) ■ (41) 

where Qk and are real functions. Then we can rewrite the formula (38) for TV = 1 in a more explicit way: 

cosai 
sinax 
sinhQi 



rs+i = r B H ^7- I sinai I . (42) 

cosher 1 ' 



The representation (41) is especially convenient in the su(2)-AKNS case. If $s = $0 (see (35)) then Qk 
and a.k are linear in x and t (x = x 1 , t = x 2 ): 

Qk = 2a;ImA fc + 2tlmw k + Qko , , . 

(43) 

ak = 2x(X — RcAfc) + 2t(u> — Rcwfe) + otko , 

where co := to(X), ujk ■= w(Afe) and Qko, afco are constant. The soliton surface ro degenerates to the straight 
line, ro = — (x + o/(A)t)e3. Therefore the formula (42), valid for any background, is especially useful in the 
case of the trivial background (35) (r B = ro, Qi and ot\ are given by (43)). In this case tb+at, denoted by 
r^v, describes the interaction of N solitons. A single soliton in the LIA case, w(A) = — 2A 2 , has been first 
found by Hasimoto [41]. 

Physical characteristics of the single soliton solution ri can be computed in the standard way. First 
of all we determine the maximum of the wave envelope (<5i = 0) which performs a helical movement. The 
group velocity vf of the soliton is computed as the velocity of the maximum along e 3 axis. The rotation 
rate of the maximum is denoted by fii. Then we determine positions of the individual wave peaks (a\ = 0). 
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Their velocity (or the phase velocity) is almost constant sufficiently far from the envelope maximum. The 
phase velocity of the single soliton wave will be denoted by v\ h . Straightforward computations yield: 

V k ~ KIT _ U W ' V k - A— RoAfc — V / ' 

(44) 

= 2(w(A) - R CWfe ) - 2(A - RcA fc )is^ . 
To describe interactions of solitons (the case TV > 1) it is convenient to introduce parameters Aj k , Sjk 

e A ]k +tS Jk ._ (Aj - A fc )(A- At) 
(Xj - Afe)(A - Afe) 

In the case TV = 2 we denote A := A12 = A21, Si := S12, 82 = S21 (compare [16]). The parameters are not 
independent. Indeed, di(e A ~ tSl — 1) = d2(e A+t52 — 1). Now we can write down a compact formula for tb+2'- 

I / di {e9 2 cosaf + e~® 2 cosa^) + d 2 (e® 1 cosaj + e~ Ql cosa. 
i"B+2 = r B + — I d\ (e® 2 sina^ + e~ Q2 sina^J + rf 2 (e Ql sinaj + e~ Ql sina." 

y d + sinlXQi + Q2) + cL sinh(<3i - Q 2 ) + d a sir^a! - a 2 ) 

where D = coshQi coshC^ coshA+sinhQi sinhQ 2 sinhA+sinhAcos(a!i — a 2 ), &f = cti±&i, d± = (di±d 2 )e ±A 
and do = di sm5\ — —d 2 sin<5 2 . The asymptotic behaviour of r 2 can be calculated easily. If 7^ v\ then we 
consider the limit Q 2 — >• ±00 (assuming | Qi | <C | Q2 !)• Thus 



2 




r 2 Q ^°°=l 1+ -g— - I sKaubJx) 



cos(ai ± <5i) 
sinTai ± <5i) 

±d 2 - z - w'(A)i / ™ MW \ sinhfQ! ± A) 

The shape and the velocity of the soliton do not change during the interaction. The only result of the 
interaction is the phase shift. In fact we have two phase shifts: the shift A p h along e 3 axis and the shift of 
the angular variable ot\. The first one is much more important and can be measured in experiments. It can 
be easily calculated: 

A 

A ph = 2d 2 + —— . 

ImAi 

The multi-soliton solutions are parameterized by complex eigenvalues Ai, . . . , Xn (in the sequel we denote 
Afe = CLk + ibk)- It is important to express the solutions by a set of 27V parameters which admit a physical 
interpretation, like dk, v g k , v% , fife, A p h etc. (compare [60]). The number of the physical parameters is much 
greater than 27V and one may expect a lot of constraints which can be checked experimentaly. The definition 
(39) suggests the following change of variables: 

. . , 1 a k ~X + ib k . . 1 Ck+idk 

c k + id k :== - = - 2 , Afe-A= — = 2 , , 2 • I 45 ) 

Afe — A (a fe - A) 2 + b% c fe - id k c{ + d% 

Thus the parameters Ck,dk are expressed by ak,bk and vice versa. In particular we have: 

e A jk+iSjk ._ (dj - d k ) + i( Cj - c k ) 
(dj + d k ) + i(cj - c k ) 

The physical meaning of dk is clear while the interpretation of Cfe is, in general, a non-trivial problem. 
However, in the case of the Localized Induction Equation (the case w(A) = — 2aA 2 ) the interpretation is 
quite simple. Namely, Cfe = v^/Clk which means that Cfe is the distance travelled by the wave envelope during 
one full circle of the wave maximum [16]. 

Let us consider the case of Localized Induction Equation with axial flow (37), w(A) = — 2aA 2 — 4/3 A 3 . We 
have w(0) = u/(0) = and the formulas (44) assume the form 

vl = 4/3(bl 3a 2 ) - aa k , vf = 2aa~ k \bl - a 2 ) + 4/3(3 2 - a 2 ) , 
Ofc = -4(a 2 +& 2 )(a + 4afc/3) . 

Using (46) and (45) we can easily check that 

2a 16/3 _ 2(v 9 k -v?) 
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Finally, let us try to parameterize the A^-soliton solution by the parameters dk and vf, which are most 
convenient from the experimental point of view [60]. It is sufficient to express Ck in terms of dk and v 9 k . By 
(46) and (45) we have: 

(4 + 4)M = Wl - 34) - ac k (ct + d\) . 

Therefore, c k is a root of the algebraic equation of the 4-th order with coefficients parameterized by dk, 
a and (3. 

Compact formulas for 7V-soliton surfaces can be derived from (38) for an arbitrary N [14]: 

i ( N ImA N N \ 

VB+N = T B + ■= [Y, | , m , fc ,2 - S S ' 

z \k=i I a - A fe | j=1 fc=1 y 

where i? := A -1 and A is NxN matrix with coefficients Ajk := i(£,k£,j + 1)(A — Aj)(A — ^fcXAj — -^fc) -1 , ffej 
are 2x2 matrices defined by 

I, W ' 

and, finally, ^ arc defined by (40). Note that the surface r^+^r is expressed solely in terms of the background 
wave function $s and 2N complex parameters: and jk '■= Pki/Pk2- Obviously, one can use Ck and dk 
instead of . 

We complete this section with few remarks on the general iV-soliton case. To compute 7V-soliton addition 
to the surface r := |a=a ^ assume the Darboux matrix in a general form 

N 

D = J2c k (\-Xo) k ■ 

k=0 

The matrices Ck are computed from the following linear system: 

N 

Y,C k {K - \o) k H^)P v -0 , (u = l,...,Nn) , 

k=0 

where A„ <G C and p„ e C™ are constant. Of course, one should take care of reductions which can result in 
some constraints on \ v , p v . The formula (21) assumes the form: f = r + < f>~ 1 C ~ 1 Ci < l> . 

Spectral problems from geometry 

It would be very important to be able to discern integrable classes of surfaces. Here we present shortly how 
to approach this problem in a natural way using Lie symmetries. Gauss- Weingarten equations, especially 
in the form (18), have very similar form to the Zakharov-Shabat linear problem (1). The only difference is 
an absence of a spectral parameter. It is well known that in many cases the spectral parameter is a group 
parameter, i.e. there exists a symmetry of GMC equations (usually a simple one, like a scaling or Lorentz 
or Galilean boost) which changes GW equations [57, 69]. The transformed GW equations contain explicitly 
the group parameter. We have developed a systematic approach to study the problem [15, 17, 24, 54]. It 
consists in computing two algebras of Lie symmetries: the algebra A of symmetries of GMC equations and 
the algebra A 1 of symmetries of GW equations. Always A' C A. Reasonable candidates for the spectral 
parameter are provided by vector fields v such that v G A and v ^ A'. For more details and references see 
[15, 21, 24], compare also [47]. 

Let us consider the following problem. We start from a given class of surfaces (i.e., GW equations are 
given). Suppose that A' ^ A. Thus using an appropriate symmetry of GMC eqs. we can insert a parameter 
into GW equations (18) to obtain a linear problem of the form (1). Then we apply the Sym-Tafel formula 
(20). 

Is the obtained class of surfaces identical with the class we started from? We have no general answer 
yet. In the following examples the answer is positive. 

• Pseudospherical surfaces. The fundamental forms (26) with £ = 1 and the symmetry u — (~ 1 u, 
v = (v of the sine-Gordon equation yield the linear problem (25). 
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• Constant mean curvature surfaces. The fundamental forms (29) with k = 0, the symmetry 
u = u cosk — v sinn and v = v cosk — u sin/t of the equation (28), and the gauge transformation 
$ = exp(«e2)$ yield the linear problem (27). 

• Spherical surfaces. The fundamental forms (31) with k — and the symmetry u = u cosk — v sinn 
and v = v cosk — u sinK of the elliptic sinh-Gordon equation (28) yield exactly the linear problem (30). 

• Bianchi surfaces. The fundamental forms (34) and the symmetry a = a/(, b = b(, f = //(l — 2A/) 
and g = g/(l + 2\g), where C = (1 - 2A/) 1 /2(i + 2Xg)-^ 2 = (1 - 2\~g) 1 ' 2 (l + 2A/)- 1 / 2 , yield the linear 
problem (32). 

In the second case it was necessary to perform a gauge transformation dependent on the spectral parameter 
(otherwise, starting from H = const surfaces, one obtains spherical surfaces) [31]. 



5 Integrable geometries and Clifford algebras 

The examples presented above are associated with the SU(2) group. Recently, the soliton surfaces approach 
has been applied in more complicated cases (although the word "surfaces" is slightly misleading: one can 
consider submanifolds of higher dimensions). It is convenient to use Clifford algebras C(p,q) generated by 
elements ei, . . . , e m (m = p + q) satisfying 

e M e„ + e v e^ = 2t]^ v , (47) 

where r]^ v = for \i ^ v, rj^^ — 1 for \i — 1, . . . ,p and r)^ = — 1 for n = p + 1, . . . , m. In the sequel the 
matrices Uj of the linear problem are linear combinations of e^e,,. Then the function 'J assumes values 
in the group Spin(p, q) which is the double covering of SO(p, q). Note that SU(2) can be identified with 
Spin(3). 

Isothermic surfaces (isothermic immersions in E 3 ) can be defined as surfaces admitting infinitesimal 
isometries preserving the mean curvature. Their curvature lines parameterized in a proper way form a 
conformal coordinate system. In other words, there exist local coordinates x x ,x 2 in which fundamental 
forms read as follows [4]: 

I = e 2i> ((dx 1 ) 2 + {dx 2 ) 2 ) , 
II = e^ikiidx 1 ) 2 + hidx 2 ) 2 ) , (48) 

where k\,k,2 and 1? depend on x l ,x 2 . Recently we found the following linear problem [20, 25] which enable 
one to study isothermic surfaces using powerfull tools of the theory of solitons: 

$,1 = i ei (—#,2 e 2 — k2e" & e 3 + A sinlr$ + A cosh?? e 5 ) $ , 

1/9 \ ( 4Q ) 

= 5 e 2 (— #,1 e 1 — fcie"e 3 + A coshi9e 4 + A sinhi? e 5 ) $ , 

where e k satisfy (47) with (17^) = diag(l, 1, 1, 1, — 1) and ie 1 e2e 3 e 4 e 5 = 1. The Sym-Tafcl formula needs 
some modification to be applied to isothermic surfaces. Namely, the formula (20) for A = defines a surface 
in 6-dim. space. Projecting the surface onto appropriate orthogonal 3-dim. subspaces we obtain a pair of 
dual isothermic surfaces. Indeed, one can prove (see [20]) that 

r± := i(lTe 4 e 5 )$- 1 $,A |a=o (50) 

are isothermic surfaces immersed in Euclidean spaces spanned by efc(e 4 ± e 5 ) (k = 1, 2, 3), respectively. The 
fundamental forms of the surface r + are given exactly by (48) while the "dual surface" r_ is the so called 
Christoffel transform of r + [4, 9]. The Darboux matrix for the linear problem (49) reads 

D = /,T o («i(Pi e i +P2e 2 +P3G3) + A(coshx e 4 + sinhx e 5 )) , 
V A + K i 

where K\ is a real parameter and X,Pi,P2 5 P3 are real functions which can be explicitly expressed by $ 
evaluated at A = [18, 20]. The corresponding Darboux-Bianchi transformation for the surfaces (50), 

r± = r± + — e ±x (±p 1 e T ' d r ± , 1 +p 2 e Ti) Y± l2 -p^±) , 
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is identical with the classical Darboux-Bianchi transformation for isothermic surfaces (compare [4]). The 
connection of the old branch of the classical differential geometry with the modern soliton theory via the 
linear problem (49) has already started a series of new interesting developments in this field [9, 12, 25]. 

Space forms of dimension n in R 2n_1 . The Lobachevsky plane can not be immersed globally in 
R 3 and only local immersions (pseudospherical surfaces in R 3 ) are possible. Analogical situation has place 
for immersions of n-space forms (spaces of constant curvature) for n > 2. There are theorems on non- 
existence of global immersions (see [33] and references cited therein) but local immersions are known (at 
least implicitly) [3, 64]. If the curvature is constant and negative then there exists a coordinate system such 
that all fundamental forms are diagonal. Moreover, the immersions turn out to be parameterized by an 
orthogonal nxn matrix function (ay) (see, for example, [1, 3, 81]). The fundamental forms read: 

n n 

I = J2 aljidxi) 2 , II m = ]T aij a mj (dxi) 2 . (51) 
i=i i=i 

where m = 2, . . . ,n and Ym=i a ji a jk = <>jk- The coefficients ay have to satisfy the Gauss-Mainardi-Codazzi- 
Ricci equations. The equations are integrable: the Backhand transformation [81], the inverse scattering 
method [1] and the loop group approach [33] have been applied. The following linear system with the 
spectral parameter A is a modification of the so(n, n) spectral problem presented by Ablowitz, Beals and 
Tenenblat ([1], see also [33, 38]): 

I ( \ 2 - I A 2 + l™ " \ 

= 2 ( 2A aijGl + 2A ^ ak ° Gk + ^ 7fe J e "+ fc ) e «+J* > ( 52 ) 

\ k=2 fc=l / 

where (ay) and (7y) are nxn matrices, -fkk = and e M satisfy (47) with ?7 M „ = 8^ v . The compatiblity 
conditions for the linear system (52) yield that (ay) is orthogonal and 

dlj,k= IkjCLlk , lij,k= liklkj , lkj,k+ljk,j + J2i=llij'Yik = dljdlk , (53) 

where indices i, j, k are distinct and, like the index I, run from 1 to n. Let us define a map F by the Sym 
formula 

F := I a=i • (54) 

It defines (at least locally) an n-dimcnsional manifold (with possible singularities) immersed, obviously, in 
the space of dimension n(2n — 1) isomorphic with the Lie algebra so(2n). One can show that in fact F 
is confined to some Euclidean space of dimension In — 1. Indeed, the unit tangent vectors are given by 
Ej = i5 ,_1 eie„ + j\E' while the normal space is spanned by Nj := ^^ 1 eiej , i> (j = 2, . . . , n). Differentiating 
Nj we obtain elements of the tangent space, Nk,j = cikjEk- The corresponding fundamental forms are given 
exactly by (51). Thus we conclude that the formula (54) is an explicit expression for the n-dimcnsional 
manifolds of constant sectional curvature locally immersed in a Euclidean space of dimension (2n — 1) [22]. 

Lame equations. We proceed to the general description of orthogonal coordinate systems in Euclidean 
spaces. The metric 

/ = Hlidx 1 ) 2 + ... + Hl(dx n f , (55) 
is induced by an orthogonal system in E ra provided that Hk satisfy the Lame equations [32] 

\~ ) ' fe + [~ ) >3 + ^=1 Hf _ U ' 



tt Hj, k H k ,j Hj, k Hj 

11 j nk— Hk T H . 



(56) 



Consider the following linear problem: 

y, j = e j (\a. j +b j )V , (57) 

where ei,e2...,e2 n generate the Clifford algebra C(2n), &j := ^ ELi a j> e n+i! D j : ~ 1 ELi ft' e > an< ^ 
f3jj = 0. The compatibility conditions for (57) imply that the matrix (ay) is orthogonal, the coefficients /3y 
are given by (3jk = —otji,k /ctki and, finally, Hk := a/y satisfy (for any fixed j) the Lame equations (56). Let 
us define 



F := y- 1 ^. 



X \x=o 
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The tangent vectors read F,j = "J 1 Uj,\ W|a=o = *o le j' a j*o, where \&o := ^(x 1 , • • • , x n ; 0) is contained in 
the subalgebra generated by ei, . . . ,e„. It is convenient to consider projections n fe defined by 

n n 

X=J2 X ii^n+i — ► U k X = X™ := J2 X jk ej , (58) 

and to apply them to F. As a result we obtain F^\F^ 2 \ . . . Differentiating F,^ = a^fc^Q 1 e J \E'o 

we derive that (F,f } | F,f ) ) = a ik ajkSij. Therefore, each map F^ has the diagonal metric tensor (55). 
Computing second derivatives we can check that F^ defines an immersion E™ — > E™, i.e., coordinates in 
E" [23] . It would be very interesting to interpret the results of [45] in our formalism. By the way, the linear 
problem (49) can also be derived from (57) as a consequence of the compatibility conditions. We just have 
to assume aj as a linear combination of e 4 , e 5 and hj as a combination of ei, e 2 , e 3 . 

6 Conclusions 

The application of soliton techniques to the differential geometry revealed deep relations between these two 
areas [7, 8, 70, 71, 79, 82]. In fact, many ideas of the soliton theory can be traced back to XIX century. I would 
like to point out that the classical differential geometry studied a lot of special immersions and interesting 
transformations between them [5, 32, 84]. It is intriguing that the Sym-Tafcl formula usually reconstructs 
immersions corresponding to a given GMC system. Then, using standard methods of constructing soliton 
solutions, one is able to derive the classical transformations of surfaces. 

The other important application of the Sym-Tafel formula is the construction of discrete surfaces. There 
is still not so clear how to define discrete analogues of integrable classes of surfaces. Surprisinly effective way 
is provided by applying the formula (20) to the corresponding discretization of the linear problem (compare 
[9, 10, 22]). " 

The recent results of Doliwa and Santini suggest even more general approach. They consider immersions 
in a sphere S n and the radius of the sphere is related to spectral parameter. In the limiting case (infinite 
radius) S n — > R™ and one can derive the Sym-Tafel formula. The results concerning integrable evolutions 
of curves are very promising [29, 30]. 
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